Ion trap transducers for quantum electromechanical oscillators 
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An enduring challenge for contemporary physics is to experimentally observe and control quantum 
behavior in macroscopic systems. We show that a single trapped atomic ion could be used to probe 
the quantum nature of a mesoscopic mechanical oscillator precooled to 4K, and furthermore, to cool 
the oscillator with high efficiency to its quantum ground state. The proposed experiment could be 
performed using currently available technology. 
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A quantum electromechanical system (QEMS) is a de- 
vice where the quantum nature of either the electronic 
or mechanical degrees of freedom becomes important in 
the observable behavior P, 0] . Potential applications of 
QEMS include single spin detection single molecule 
mass spectrometry 4], and readout for quantum infor- 
mation devices 5] . An excellent example of such a sys- 
tem is a rf single electron-transistor integrated with a 
low loss, high frequency nanomechanical resonator which 
demonstrated both continuous position detection approx- 
imately a factor of five away from the uncertainty princi- 
ple limit, and direct observations of the mechanical mode 
occupation factor as low as NwSS However, this de- 
vice suffers from both practical and fundamental limi- 
tations; picowatt levels of dissipation at the transistor 
is suspected to have blocked passive cooling of the me- 
chanics to lower temperatures and occupation factors, 
as well as the fundamental sensitivity limitation due to 
coupling to an intrinsically non QND (quantum non- 
demolition) variable such as position. Other methods 
to perform measurements on quantum limited mechan- 
ical systems are clearly needed. Furthermore, coupling 
atomic systems to nano-electronic and nano-mechanical 
devices appears to be a very promising and exciting new 
frontier, where it is hoped one might combine the un- 
matched quantum coherence and detection efRciencies of 
atomic physics, with the sophisticated electronic and mi- 
crostructures which are possible in the condensed matter 
realm 

In this paper we investigate the use of a laser-cooled, 
trapped atomic ion to both monitor and manipulate the 
number-state of a QEMS oscillator, following earlier sug- 
gestions by Wineland and others 0, 0|. As the tem- 
perature of an ion's vibrational degree of freedom can 
be < 10"'^ K, it is hoped that such a transducer could 
generate far less thermal power than low temperature 
single-electron devices. Furthermore, due to the excellent 
optical readout achieved in ion trap systems via fluores- 
cence shelving [ill ], the ion trap transducer is expected 



to posses sensitivity of QEMS energy at the level of a 
single quanta, a detection which is very difhcult with a 
simple linear coupling to displacement. Finally, we show 
that the ion-QEMS system can be configured to provide 
a very effective cooling mechanism for the mechanical 
oscillators Q and estimate that it should be possible to 
cool a QEMS cantilever, precooled only to 4 Kelvin, to 
its ground state with very high efficiency. 

Here we describe a quantum model of a single trapped 
atomic ion which is electrostatically coupled to a very 
small doubly clamped cantilever. Fig. ^ This coupling 
can be switched on and off using an external bias voltage 
at an electrode on the oscillator. The ion is held in a 
mesoscopic micro- fabricated quadrupole Paul trap O, 
and laser cooled using resolved sideband cooling ISj. 
External lasers are used to couple the internal electronic 
states of the ion to its vibrational degree of freedom. The 



FIG. 1: Schematic representation of a single trapped atomic 
ion coupled to an nanoelectromecanical oscillator. 

bias gate on the oscillator carries a charge Q = CoVo{t) 
where Co is the capacitance of the gate. We allow for 
the possibility for the bias gate voltage Vo{t) to be time 
dependent so that it may be set to zero to turn off the 
electrostatic coupling to the trapped ion. The ion carries 
charge +e. In the geometry of Fig. ^ the interaction 
energy between the ion and the oscillator is given by 
Vc = keVoCo/\d + X{t) — x{t)\ where d is the equihbrium 
separation of the oscillator centre of mass position and 
the ion. We set the equilibrium position of the ion at the 
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origin for simplicity. In that case X{t) and x{t) represent 
the small oscillations of the QEMS oscillator and the 
ion around their equilibrium positions, respectively. 
We now assume that the deviations from equilibrium 
are small compared to the equilibrium separation 
d and expand to second order in {X{t) — x{t))/d. 
The interaction energy is then given by Vc = 

{keVoCo/d) \l - {X{t) - x{t))/d + {X{t) - x{t)f/d^ . 

The linear term may be absorbed into the definition of 
the equilibrium positions. The quadratic term includes 
a coupling of the oscillators and a renormalisation 
of the oscillation frequency for both the ion and the 
oscillator. This frequency shift is a small perturbation 
of the bare frequencies and is neglected. We then see 
that the interaction Hamiltonian coupling the ion to the 
oscillator is = —x^{t)x{t) where x — "^^k eVJ Q n I <P ■ 
A similar ion-oscillator system is described in |14j . 

In ion trap technology a laser can be directed onto the 
ion to induce transitions between the internal electronic 
states (ij. The electric field seen by the ion depends 
on its motion in the trap, thus the external laser can 
couple the internal electronic state to the vibrational de- 
gree of freedom. In this way information on the vibra- 
tional motion can be transferred to the electronic degree 
of freedom where it can be efficientl y rn easured using the 
technique of fiuorescence shelving [llj. Such measure- 
ments are very nearly perfect projective measurements 
onto one of two electronic states, {\g), |e)}. If the laser 
is detuned to the first red (blue) motional sideband, the 
ion can make a transition from the ground to the excited 
state absorbing a laser photon and a single quanta of vi- 
brational en ergy is simultaneously added (deducted) in 
the process [iSjj. The probability of excitation depends 
on the phonon distribution in the trap, so subsequent 
readout of the electronic state then reveals information 
on the phonon number in the trap. 

The overall scheme for using the trapped ion as a trans- 
ducer of QEMS motion is as follows. The ion is cooled to 
the vibrational ground state and prepared in some appro- 
priate electronic state. Next, a voltage is applied to the 
oscillator gate, coupling the motion of the QEMS to the 
vibrational state of the ion. This transfers phonons from 
the oscillator to the ion. In the next step an external 
laser couples the trap phonons onto the electronic de- 
gree of freedom which is then readout using fluorescence 
shelving. 

We define dimcnsionless annihilation and creation op- 
erators for the ionic vibration and the QEMS vibration 
using the ground state standard deviation of position of 
the oscillator as a convenient length scale in each case, 
X = {hl{2Muj)f'^ {a+a^) and x = {h/{2mv)f''^ {b+b^) 
where v is the vibrational frequency of the ion in the trap 
(also called the secular frequency), u is the resonant fre- 
quency of the QEMS oscillator, and M, m are the masses 
of the QEMS and the ion respectively. The total Hamil- 



tonian for the three coupled systems may be written as 
H = huja^a+hvb%-TiK{a+a^)\b+b'^)+h^<7-,+He where 
describes the interaction between the external laser 
and the electronic states of the ion (see below) . The cou- 
pling constant k is given 

The electronic transition frequency is oja and Uz — 
|e)(e| — \g){g\ is a Pauli matrix. The interaction between 
the control laser and the electronic transition is given in 
the dipolc approximation by He = 2Mlax sm.{kix — uJit) 
where LOi and ki are the frequency and wave vector of the 
laser, J7 is the effective Rabi frequency for the transition 
and da; = |e) (g| -|- \g) (e| . We will work in the Lamb-Dicke 
limit in which the amplitude of the ion's motion in the 
direction of radiation is much less than the wave length 
of the laser yjub + Ir/ = y/nb + Ih {h/{2miy))^^^ « 1 in 
which case we may take = 2hr]n(Tx{b + 6^) sm{uJit). 
We now move to the interaction picture and make the 
rotating wave approximation for both the QEMS-ion 
coupling and the electronic coupling, and assume that 
uji — LUA — h' or llH = LOA + I' so that we are resolving 
the first red or blue sideband. The total Hamiltonian is 
Hj — hAa^a — hn^a^b + a6^) + iJgf, with A. = uj — v 
and Hsb — ^.9(&o'+ + &^o'+) for the red sideband and 
Hsb — figih^ a j^+ba\_) for the blue sideband where g = rjQ, 
and CT+ — |e) (^j. In our model both k and g can be turned 
on and off. The measurement protocol has two stages. 
In stage I, the oscillator and cantilever are coupled for a 
time T by setting k ^ 0, g = 0. In stage II, k, is turned 
off and the electronic and vibrational degrees of freedom 
of the ion are coupled by red and blue sideband excita- 
tion. At the end of this stage, the ionic electronic state 
is measured. 

Using the ion as an ultrasensitive sensor for the can- 
tilever motion is experimentally realistic. Assuming a 
19.7 MHz cantilever, as reported recently by LaHaye et 
al. _;6;j, a radio-frequency microfabricated ion trap [l^ 
could be used to confine the ion at a distance (3 of 50/ito 
from the cantilever at a secular frequency of 19.7 MHz. 
Applying a static voltage of 7.5V one would obtain a 
coupling frequency 27r x 52.5 kHz for a cadmium 

ion. Assuming the cantilever at a temperature at 4 K 
(liquid Helium), the cantilever contains on the order of 
1000 quanta of motion. Therefore one would only re- 
quire an interaction time of 5 /xs for the ion to undergo, 
on average, a single phonon excitation. Anharmonicity in 
the cantilever motion is included in the measured Q (too 
small to observe at 4K) and will not measurably distort 
the cantilever lineshape for any reasonable temperature 
(isf . An upper bound for the width of the secular mo- 
tion frequency due to the anharmonicity of the ion mo- 
tion at fib = 4000 is found to be Az/ « v{z/ fiy^ ~ where 
z = {fifib/ {2mv)Y/'^ is the extent of the atomic wave func- 
tion and (3 is the dimension of the ion trap. For the 
parameters considered here, Az/ <C k, so ion trap anhar- 
monicities are not expected to degrade the ion-cantilever 
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coupling. The QEMS oscillator is always coupled to a 
thermalising reservoir and before the coupling to the ion 
is turned on we assume that the oscillator is in a ther- 
mal equilibrium with mean vibrational quantum number 
rxao = (e''"/'=«^ - ^ kBT/ihu;). 

When the coupling between the oscillator and the ion 
is turned on, the ion becomes indirectly coupled to the 
phonon reservoir of the QEMS oscillator. This may be 
modelled using a quantum optics master equation for the 
oscillator (if the frequency of the oscillators is not too 
low) that describes the joint density matrix operator R{t) 
for the QEMS oscillator and the ionic vibration motion 
during stage I, 

dR 

= ~iA[a'^a,R]+iK[a''b + ab^R] 

at 

+7a {nao + 1) 'D[a]R + -ianaoV[a)]R 

+HlV[b]R + tJL2V[b^]R, (1) 

where V[0]p = OpO"^ - {0'^0p + pO'^0)/2 is defined for 
arbitrary operators O and p, and the damping rate 7^ 
is related to the quality factor Q of the QEMS oscilla- 
tor by the expression 7q ~ lo/Q. Note that choosing the 
damping rates /ii = 7b(nbo + 1), ^J'■2 — Jbfibo corresponds 
to the dynamics for the ionic vibrational motion induced 
by interacting with a thermal bath, as in the case for the 
QEMS oscillator. However, if the main noise source cou- 
pled directly to the vibrational modes of the trapped ion 
is the fluctuating electric field generated in the trap elec- 
trodes , we can model it as a classical stochastic elec- 
tric field. In this case, we take pi = u-?. = (1/ti), where 
Ti is the characteristic heating time [13 • Extrapolated 
from recent experiments and associated theory for char- 
acteristic heating times we obtain an expected 
approximate heating rate on the order of 0.06 quanta per 
ms for the above geometry for the cadmium ion. Compar- 
ing this heating rate with the coupling constant k from 
our example above it is reasonable to neglect ion heating 
during the cantilever-ion interaction. We are interested 
in the mean phonon number nh{T) — (6^(r)6(T)) of the 
thermal state of the ionic vibrational modes at the end 
of stage one. Using Eq. the equations of motion for 
quadratic moments form a closed set and do not involve 
higher-order moments. 

Neglecting heating due to a classical stochastic field 
during the ion-cantilever interaction (as justified above) 
and assuming the secular ion frequency equal to the reso- 
nance frequency of the cantilever one can derive a simple 
analytical expression {A — pi — p2 = 0): 

(n^o - n-ho) 



«6(t) 



naO 
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(2) 



stage one the QEMS oscillator and ionic vibrational 
mode are in thermal states with mean phonon num- 
ber fiaQ and nfto, respectively. For 7a <C k one ob- 
tains nfc(T) — fiao — {fiao ~ nf,o)e~'^'''^^^cos^{KT). The 
mean oscillator phonon number can also be found to be 
riair) = Hao - {fiaO ~ "-f,o) e"''"'^/^ sin^ {^■yT) /^'^. In 
Fig. 21 we plot nb(r) (obtained from Eq. 101) using re- 
alistic parameters corresponding to the mechanical can- 
tilever of Ref. 0. Assuming fiao = 4000, fihQ = 0, A = 0, 
K = 27r X 52.5 kHz and assuming a Q factor for the can- 
tilever of 30000 we find that the effect of the interaction 
of the cantilever with the thermal bath is approximately 
negligible during the first quarter coupling period, how- 
ever, becomes significant on longer timescales. In the 
limit << 1 we find that n.(,(r) « fiaQ s\t?{kt). 



n(T): 




where Vl^ = -^/k^ — (70/4)^. The initial conditions 
used to obtain Eq. ||2Jl assume that at the start of 
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FIG. 2: The mean phonon number of the ion fibij) as a func- 
tion of time r obtained from Eq. ^ is plotted (solid line) 
and compared to the case where the couphng of cantilever 
and thermal reservoir is neglected during the cantilever-ion 
interaction (dashed line). The mean phonon number of the 
cantilever fiaij) is plotted as dotted line. The inset shows 
quadratic short time behaviour of the ion. 

From Fig. |21 we see that eventually the ionic vibra- 
tional degree of freedom comes into thermal equilibrium 
with the cantilever. However due to the very strong cou- 
pling between the two, the dynamics of this process is 
certainly not exponential. The oscillations in Fig. [3 
are thus evidence of non-exponential relaxation. The 
experiment proposed here provides a convenient way to 
study the transition from exponential relaxation to non- 
exponential short time behavior by varying the coupling 
strength k with respect to the damping rate of the can- 
tilever, "fa- 

The readout process of stage II is similar to the exper- 
iments described in Ref. 0]. In stage II we couple the 
electronic state of the ion to its vibrational motion for a 
time T using the first red and blue sideband transitions. 
If we write the probability for the atom to found in the 
excited state after time T as P^{T) and P^iT) for red 
and blue sideband excitation respectively it can be shown 
0,01 that the mean phonon number ^^(t) is given by 
nb(r)/(l + fi6(r)) = P^{T)/P^{T) = R,. As this re- 
sult is independent of the coupling time T, noise due to 
heating during stage two can be kept small by minimiz- 
ing T. Indeed such experiments are routinely performed 
|l8l nil l2oj and a phonon number smaller than unity 
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can be determined. As the parameters k, r are assumed 
to be known this gives fiao directly. In particular, for 
the short time region of Fig. |21 (and << 1) we find 
that Pf(T)/P^{T) « Uaon'^T'^- Thus measurement of 
the ratio of excitation probability on the first red and 
blue sideband yields fiao directly. Note that rib(r) should 
be on the order or smaller than Umax — 20 for reliable 
measurement fiolll^lll^ . 

From Fig. El we see that it is possible to almost com- 
pletely exchange the mean phonon number of the can- 
tilever and the ionic vibration. Since the ion begins in the 
vibrational ground state, the result after one exchange 
time is to transfer the thermal energy of the cantilever 
rapidly to the ionic vibration, leaving the cantilever near 
its ground state. 

A number of cooling procedures appear possible. One 
could simply dump the resulting hot ion or one could im- 
plement two ion traps that are located adjacent to the 
cantilever: one ion trap to be used to cool the mechan- 
ics to the ground state, while the other ion trap is used 
to perform the quantum measurement of the cantilever. 
Alternately, an iterative cooling mechanism could be ac- 
complished by decoupling the resulting ion from the cold 
cantilever (instead of dumping it) by detuning the ion 
trap secular frequency and then laser cooling the ion. 
One could also accomplish cooling of the cantilever and 
obtain a reasonable temperature estimate using just a 
single ion provided fiaoT^la/i^ < nmax (meaning that the 
ion carries less than Umax phonons after one full cou- 
pling period). In this scheme one waits for a full cou- 
pling period to map out the minimum of ^^(t). Using 
this revival one could deduce the optimum cooling in- 
teraction time (half the coupling period) and the asso- 
ciated temperature of the cantilever. Although the cold 
state is metastable for all of the schemes above, with 
a lifetime determined by 7a, we expect a significantly 
lower resonator temperature compared to what has been 
achieved with continuous refrigeration of the cantilever 
bath, fiao ~ 50 0. Finally, one could couple the ion 
continuously resulting in sympathetic cooling and ob- 
tain a stable low temperature of the cantilever. This 
can achieve a mean phonon number of the cantilever less 
than unity [l^ . 

Approaching the mechanical ground state may be facil- 
itated by starting with a much higher frequency resonator 
(eg. a 1 Ghz cantilever ^) and dilution refrigeration. 
While it may be difficult to fabricate an ion trap with 
such high secular frequency, it is possible to couple the 
cantilever to the micromotion [lol |23 | in the ion trap 
that occurs at Q^f + and Qrf — v where 11^/ is the 
radio-frequency driving frequency of the ponderomotive 
ion trapping potential. Micromotion could then be cou- 
pled via laser excitation j2^J into the electronic state of 
the ion. 



This scheme is not a single shot readout of phonon 
number in the oscillator, rather it enables a statistical in- 
ference of the mean phonon number of the oscillator. Can 
it be used to measure a weak classical force ? The sensi- 
tivity to changes in fihij) is best for values of fif,(r) less 
the unity. This can be achieved by ensuring fiaoK^T^ < 1. 
It should then be possible to detect a change in the in- 
ferred mean phonon number of the oscillator of the order 
of one quanta. If a weak classical force is continuously 
applied to the cantilever after it has been prepared close 
to the ground state (eg. via sympathetic cooling via the 
ion) , its equilibrium phonon distribution will shift by an 
amount proportional to the square of the ratio of applied 
force to the energy damping rate. Thus we should be 
able to infer the size of a classical force so weak as to 
shift the mean phonon number of the oscillator by one 
quanta. This corresponds to a displacement sensitivity 
at the standard quantum limit Axsql = (?i/2mLj)^/^. 
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